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I.  INTRODUCTION 


The  behavior  of  a  composite  laminate  depends  on  variety  of 
characteristics  including  stiffness,  strength  and  behavior  under 
environmental  changes.  The  large  number  of  parameters  and  the 
extensive  amount  of  calculations  involved  in  the  characterization  of 
composite  laminates  suggest  the  use  of  electronic  computers. 

Algorithms  for  solutions  of  laminate  problems  in  various  computing 
facilities  are  given  in  [1-3]*.  In  the  present  work  an  algorithm  for 
the  solution  of  the  general  laminate  shown  in  Figure  1  is  proyided 
using  an  Apple  computer. 

A  review  of  relevant  equations  is  provided  in  Section  II  which 
includes  modulus  and  compliance  analysis,  hygrothermal  effects,  strain 
computation,  and  strength  analysis.  This  material  is  based  on  a  book 
by  S.  W.  Tsai  and  H.  T.  Hahn,  [4]. 

Instruction  for  program  running  and  control  is  given  in  Section  III. 
This  includes  data  input  procedure  and  printout  control. 


*  1.  S.  W.  Tsai,  R.  Aoki,  "TI-59  Magnetic  Card  Calculator  Solutions 
to  Composite  Materials  Formulas",  AFML-TR-79-4040. 

2.  Som  R.  Soni,  "A  Digital  Algorithm  for  Composite  Laminate 
Analysis-Fortran" ,  AFWAL-TR-81-4073. 

3.  Won  J.  Park,  "Radio  Shack  TRS-80  Pocket  Computer  Solutions 
to  Composite  Materials  Formulas",  AFWAL-TR-81-4074. 

4.  S.  W.  Tsai,  H.  T.  Hahn,  "Introduction  to  Composite  Materials", 
Technomic  Publishing  Co.,  Westport,  CT  06880,  July  1980. 
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II.  REVIEW  OF  EQUATIONS 


A  short  review  of  relevant  equations  is  given  in  this  section. 
For  a  detailed  derivation  the  reader  is  referred  to  reference  4. 

1.  Modulus  and  Compliance  Analysis 

With  deformation  prescribed,  the  effective  loads  are  found 

from* 


N1 
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where  and  ki ,  i  =  1-3,  are  mid-surface  strain  and  curvature 
components,  respectively,  and  N^,  ,  i  =  1-3,  are  average  forces  per 

unit  length  and  average  moments  per  unit  length,  respectively.  Referring 
to  Figure  1  for  notation,  the  stiffness  matrices  in  (1.1)  are  given 
by 
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*  For  convenience  the  axis  "6"  in  reference  4  is  replaced  in  this  work  by  "3" 
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where 
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Qxx  =  moEX’  Ex  =  lon9ltuc,inal  Young's  modulus 
Qyy  =  mQEy,  Ey  =  transverse  Young's  modulus 

Qxy  =  Qyx  =  moEyvx5  vx  =  lon9itudinal  Poisson's  ratio  (1.5) 

Qss  =  Es,  E$  =  longitudinal  shear  modulus 

"0  ■  '/('  - 

With  the  aid  of  (1.1),  the  deformation  can  be  expressed  in  terms 
of  effective  loads 


where  g  =  A'1  -  gBA"1 ,  §  =  -A_1B6,  6  =  (D  -  BA^B)"1  (1.7) 

It  is  possible  to  normalize  (1.1)  and  (1.6)  with  respect  to  the  total 
laminate  thickness,  h.  The  results  are: 
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(1.10) 


2.  Hygrothermal  Analysis 

The  effective  loads  generated  by  temperature  change,  AT,  and 
moisture  content  change,  C,  are  determined  using  the  following 
procedure: 


(i)  The  nonmechanical  strain  components,  e.,  are  given  by 

ei  =  o^AT  +  B^C,  i  =  x,y  ,  es  =  0  (2.1) 

where  and  Bi  are  coefficients  of  thermal  expansion  and 
swelling,  respectively. 

N 

(ii)  The  stresses  required  to  produce  these  strains,  o-,  are 

j 

found  from 

=  Qjkek’  J,k  =  x’y  »  °s  =  0  (2.2) 

where  the  superscript  "N"  has  been  assigned  to  indicate 
nonmechanical  stresses. 
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(iii)  The  on-axis  stresses  in  (2.2)  can  be  transformed  to 
off-axis  stresses  using  (2.3) 


(iv) 


N 

N 

P 

q 

N 

°1 

1 

cos26 

N 

°2 

1 

-cos20 

N 

°3 

sin20 

(2.3) 


PN  -  (o*  ♦  o*)/2.  q" 


(cJJ  -  oJ)/2 


where 

The  effective  nonmechanical  forces  and  moments  are 
given  by 

h/2 

[N-,  M}]  =  /  4>  •  o-  [l,z]  dz,  i  =  1  -  3 

-h/2 


(2.4) 


3.  Strain  Analysis 

The  object  here  is  to  determine  on-axis  and  off-axis  interlaminar 
strains  from  prescribed  loadings  (mechanical  or  nonmechanical). 

Assuming  a  linear  strain  variation  across  the  laminate  thickness, 

i  .e. 


e  =  e  +  zk 


(3.1) 
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and  using  (1.6)  in  (3.1),  the  off-axis  strains  at  z  is  given  by 


£  *  &N  +  §  M  +  z  (g  TN  +  5  M)  (3.2) 

Next,  the  on-axis  strains  are  found  using  the  transformation  in 
(3.3) 
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-cos28 
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2cos28 

where  p  =  (e]  +  e2 )/2 ,  q  =  (e1  -  e2)/2’  r  =  C3^Z 
4.  Strength  Analysis 

In  this  work  laminate  strength  is  examined  using  two  failure 
criteria  ,  i.e.  the  Tsai-Hill  and  the  Maximum  Strain. 

In  the  maximum  strain  criterion  failure  is  assumed  when  one  of  the 
six  conditions  below  met  first 


Cy,  e$>0.  (ex,  £y,  e$)  i  anowecj  "  (X/E> 


V/Ey,  S/Es) 


(4.1) 

V  Ey,  Es  <  0:  -  "  -  =  (-X'/Ex,  -Y'/Ey,  -S/E,) 

where  X  and  X'  are  longitudinal  tensile  and  compressive  strength, 
respecti vely,  Y  and  Y'  are  transverse  tensile  and  compressive 
strength,  respectively,  and  S  is  the  shear  strength. 


Defining  strength  ratio  R  as 


^  “  £i  | all owed^l  |  imposed’  ^  "  x»  5  ^4>2) 

and  assuming  nonmechanical  strain  as  well  as  mechanical  strain  exist, 
then,  with  superscript  "M"  assigned  for  mechanical  strain 

Ei [allowed  =  R  £i  +  £1  -  eT  <4'3> 

using  (4.3)  in  (4.1),  one  has 


min. 


/  X  N  ,  Q  ,  ,  M  / Y  N  .  \  /  M  ,S  Nw  M 

(  £-  ‘  £x  +  eX)/£X  *  '£ - £y  +  ey)/Cy  *  (f - £S^eS 

x  y  s 


(4.4) 


y  y"_y 

where  X,  ?,  5  =  X,  Y,  S  for  positive  ,  i  =  x,  y,  s 
and  X,  ?,  $  =  -X1,  -Y1,  -S  for  negative  ,  i  =  x,  y,  s 
In  the  Tsai -Hi  11  criterion  failure  occurs  when 

Gij  ei|a11owed  ej|allowed  +  Gi  £i]a11owed  x>  s  ^-5) 

where  the  nonvanishing  terms  in  (4.5)  are 


Gi  =  Fj  Qij 

Gkl  =  Fij  Qik  Qjl’  1’  k>  1  =  x*  y  (4-6) 

Gss  *  "Vs>2 

Fx  =  1/X  -  1/X\  Fy  =  1/Y  -  1  / Y 1 

Fxx  *  V(XX').  Fyy  -  1/(YV),  Fxy  •  FJy  (FxxFyy)1/2  (4.7) 

Introducing  (4.3)  in  (4.5),  one  finds  two  roots  for  R,  one  positive 
and  the  other  negative.  Only  the  positive  solution  is  given  (the 
negative  root  corresponds  tc  a  reverse  straining). 
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Ill  PROGRAM  CONTROL 


The  program  language  is  in  "Applesoft"  ("BASIC"  with  some  additions) 
and  it  is  described  in  the  Apple  instruction  manual.  The  program  flow 
diagram  is  shown  in  Table  I.  Terminology  for  inrut  and  ouput  data 
is  given  in  Table  II  and  computer  memory  allocation  in  Table  III. 

The  program  listing  and  illustrative  examples  are  shown  in  page  15 
and  21,  respectively.  Program  control  and  data  input  procedure  are 
summarized  below. 

1.  Running  the  program 

With  the  disk  inserted  into  the  disk  drive,  the  program  "composite" 
is  loaded  automatically  into  the  computer  memory  once  the  computer  is  turned 
on.  Note  that  the  disk  contains  a  subprogram  used  for  printing  data  in 
scientific  format.  This  program  is  also  automatically  loaded  into  the  computer 
memory . 

2.  Data  Input  Procedure 

Data  are  inputted  through  both  program  line  editing  (before  running 
the  program)  and  computer  keyboard  during  program  run,  according  to 
the  procedure  outlined  in  Table  I. 


For  convenience,  material  properties  for  five  composites  and 
aluminum  are  stored  in  the  program  according  to  the  following  scheme: 


Program  Line 

Material  Type 

Material  Identification 

40 

40 

T300/5208  (graphite/epoxy) 

50 

50 

B(4)/5505  (boron/epoxy) 

60 

60 

AS/3501  (graphite/epoxy) 

70 

70 

Scotchply  1002'  (glass/epoxy) 

80 

80 

Kevlar  49/epoxy  (aramid/epoxy) 

90 

90 

Aluminum 

Material  selection  is  achieved  through  keyboard  by  inputting 
the  material  type  number  in  the  table  above.  Other  materials  can  be 
analyzed  by  Introducing  appropriate  material  properties  in  either  program 
line  40  to  90. 

Mechanical  forces  and  moments  on  a  per  unit  length  basis  are 

inputted  in  program  lines  940  and  950,  respectively.  The  current  values 

are  N-j  =  1/1 0^  M-GPa,  N2  =  N3  =  M1  =  M2  =  M3  =  0.  The  strength  parameter 

F*  is  inputted  in  line  1580.  Its  current  value  is  -0.5. 
xy 

3.  Printout  Control 

If  a  "hard  copy"  printout  is  desired  the  printer  should  be  activated 
prior  to  running  the  program.  The  display  and  printing  format  requires  that  both 
the  CRT  screen  and  the  printer  page  width  should  be  set  to  at  least  80  character. 

For  some  applications  a  printout  of  all  output  data  blocks  indicated 
in  Table  I  may  be  excessive.  A  selective  output  printout  is  possible 
using  the  CN(I),  I  =  1-8,  array  in  program  line  20,  as  described  in 
Table  I.  For  instance,  if  in  program  line  20  we  have  "CN(1)=0:  CN(2)=0: 

CN(3)=0:  CN ( 4 ) =0 :  CN(5)=1:  CN(6)=0:  CN(7)=0:  CN(8)=0",  then  only 
on-axis  strains  will  be  printed.  The  current  values  are  CN ( I ) =1 , 

I  =  1-8. 

4.  Program  Pause 

As  indicated  in  Table  I,  after  each  block  printout  the  program 
pauses  to  give  the  user  an  ample  time  to  observe  the  output  on  the 
CRT  screen.  This  is  accompanied  by  cursor  flashing.  To  resume  computer 
operation  press  the  "RETURN"  key.  To  eliminate  program  pause  the 
user  should  delete  the  "Get  G$"  statements  in  program  lines  2130  and 
2190. 
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-  ON-AXIS  COORDINATES 

-  LAMINATE  (OFF-AXIS)  CO. 


2 


Z 


TABLE  I  -  FLOW  DIAGRAM 

notation 


PROGRAM  LINE  DATA  INPUT 


-  KEY  BOARD  DATA  INPUT  (TYPE  THE  INPUT  VALUE  AND  THEN  "RE TURN") 


/TTTTTN  ,  .  8  -  DATA  PRINTOUT  CONTROL  PARAMETERS  GIVEN  IN 

1-1-8.  PROGRAM  LINE  20.  CN<I)=1  FOR  PRINT,  CNlII-O  FOR  NO  PRINT 

7 1  -  DATA  BLOCK  PRINTOUT 
Ik _ £ 

<0_  -  PROGRAM  PAUSES.  PRESS  "RETURN"  TO  CONTINUE  /f  f 

f  y' 

NOTE.  SEE  FIGURE  1  AND  TABLE  H  FOR  OTHER  NOTATION  (  ax'°y 


PROGRAM  LINE  40-90, 


Ex.Ey.Es.P, 

a  x,Oy.01(,0y 

X.x'.y.y'.S 


TABLE  II 


TERMINOLOGY  FOR  DATA  INPUT  AND  OUTPUT 
Data  Input 

C  (C)*  -  moisture  content 

h  ,  hc  (HO,  HC)  -  ply  and  core  thickness,  respectively,  (M)** 

N,  M  (N,  M)  -  number  of  plies  in  laminate  and  number  of  plies  below 
core,  respectively 

N^,  M^j1  (N(I,C),  M(I,0),  I  =  1-3)-  effective  mechanical  force  and 

effective  mechanical  moment  components 
(on  a  per, unit  length  basis),  respectively, 
(M-Gpa,  M-Gpa) 

Material  Type  (MT)  -  see  table  in  page  8 

E  ,  E  (EX,  EY)  -  longitudinal  and  transverse  Young's  modulus,  respectively, 
y  (Gpa) 

Es  (ES)  -  longitudinal  shear  modulus,  (Gpa) 

5  (S)  -  longitudinal  shear  strength,  (Gpa) 

X,  X'  (X,  XC)  -  longitudinal  tensile  and  compressive  strength,  respectively, 

(Gpa) 

Y,  Y1  (Y,  YC)  -  transverse  tensile  and  compressive  strength,  respectively, 

(Gpa) 

F  *  (FXY  STAR)  -  Parameter  related  to  material  strength  (See  (4.7)) 
xy 

a  ,  a  (AX,  AY)  -  coefficient  of  thermal  expansion  along  x  and  y  direction, 
y  respectively,  (1/K°) 

6  ,  8  (BX,  BY)  -  swelling  coefficient  in  x  and  y  direction,  respectively 

x  y 

v  (PX)  -  longitudinal  Poisson's  ratio  =  -c  /c 
X  y  x 

AT  (DT)  -  temperature  difference,  (K°) 

(0(1))  -  orientation  of  ith  ply  (Deg.) 

Data  Output 

A,  B,  D  (A,  B,  D)  -  stiffness  matrices,  (M-Gpa,  M^-Gpa,  M3-Gpa) 

A*,  B*,  D*(A*,  B*,  D*)  -  normalized  stiffness  matrices,  (Gpa) 

♦Quantities  in  parenthesis  indicate  program  variables 

♦♦Dimension  :  M  -  Meter,  pa  -  paschal! ,  Gpa  -  109pa,  K°  -  deg.  Kelvin 
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TABLE  II 

TERMINOLOGY  FOR  DATA  INPUT  AND  OUTPUT  (CONTINUED) 

M  fci 

Nn(N),  Nn(NN)  -  Mechanical  and  nonmechanical  force/unit  length  Vector, 
respectively,  (M-Gpa) 

mH(M),  mn(mn)  -  Mechanical  and  nonmechanical  moment/unit  length  Vector, 
respectively,  (M  -Gpa) 

N*^(N*),  N*N(NN*)  -  normalized  mechanical  and  nonmechanical  force/unit 

length  Vector,  respectively,  (Gpa) 

M  N 

M*  (M* ; ,  M*W(MN*)  -  normalized  mechanical  and  nonmechanical  moment/unit 

length  Vector,  respectively,  (Gpa) 

a,  j3,  gT,  6  (ALPHA,  BETA,  TRBETA,  DELTA)  -  compliance  matrices, 

(1/M-Gpa,  1/M2-Gpa,  1/M2-Gpa,  1/M3-Gpa) 

a*,  £*>  §*T,  (ALPHA*,  BETA*,  TRBETA*,  DELTA*)  -  normalized  compliance 

matrices,  (1 /Gpa ) 

R  (R)  -  strength  ratio 
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TABLE  III 


COMPUTER  MEMORY  STORAGE 

Modulus  and  Compliance  Components 

X(I,  J,  0)  «4— ►  A*(I,  J,)  I,  J  -  1-3 
X ( I ,  J,  1)  ►  B* 

X(I,  J,  2)  «<-►  D* 

X(I,  J,  3)  ►  a* 

x(i,  j,  4)  ►  §* 


x(i 

,  J, 

5)  «4 

-►rT 

X(I 

,  J, 

6)  -4 

-*■  «* 

Strain  Components 

(k. 

I,  o) 

- ► 

M,T^ 
e  (I) 

in  the 

kth  ply 

(k, 

I,  1) 

-4— 

- ► 

£  (i , 

in  the 

kth  Ply 

I  =  1-3  :  on-axis  strain  components  at  lower  ply  surface 
I  =  4-6  :  on-axis  strain  components  at  upper  ply  surface 
I  =  7-9  :  off-axis  strain  components  at  lower  ply  surface 
I  =  10-12  :  off-axis  strain  components  at  upper  ply  surface 

Strength  Ratio 

Rl  (k,  o)  kth  ply,  Rjsa -j  _H-j  1 1  •  lower  surface 

Rl  (k,  1)  kth  ply,  ^Tsai-Hill ’  upPer  surface 

RM  (k,  o)  kth  ply,  RMdx  strdin'lower  surface 

RM  (k,  1)  kth  ply,  RMdx  strdin, upper  surface 
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IV.  PROGRAM  LISTING 


1 :  I  N  (6> 


i :  on  •  7 ) 


1  :  ' 


1 0  D I M  CM ( I 0 ) 

20  CN  ( 1 )  *  1 :  CN  (2)  «  !:CN<3)  *  1:CN(4)  *  1;CN(5> 

N  ( 8 )  =  .1 

30  INPUT  "MATERIAL  TYPE,  MT«";M7 

10  EX  181 5  EY  -  1  3:  PX  »  .29:ES  =  7.  17:  NO  ■  .  0001 25:  X  «  1.5:  XC  =  1.5:  / 

.04:  VC  ~  .246:8  -  .068:  AX  -  .02  /  IE6:AY  =  22.5  /  1E6:BX  =»  0:  BV  = 

IF  Ml  -•=  40  THEM  GOTO  100 
5  0  EX  ■  204 5  EY  ”  18. 5:  P,\  —  .  ?3:ESr 

.061 :  YC  -  -  702:6  -  .067: AX  -  6.1  /  1E6: AY 
O:  IE  MT  =  50  THEM  GOTO  100 
00  EX  ■  138:  EY  =  9 . 96 : \  --  .  3:ES  *  7.1:  HO  ~  .000125:  X  -  1.447:  KC  ~  1.  447:  y 

-  .051 7::  VC  ”  .  2c»6:S  =  .09  5:  AX  ~  --  .3  /  1 E6 :  AY  -  28.1  /  1E6:  BX  =  F 

Y  .44:  IF  Ml  -  60  THEN  GOTO  100 
EX  -  38. 6:EY  =  8. 2  7:PX  =  .  26:  ES  ■-  4.14:  HO  •=  .  0001  25:  X  -  1 . 02.*j :  > 

.031:  YC  *  .113:5  =  .072:  AX  =  8. 


1.26: XC  «  2.5: 
30.2:  /  1E6:BX  =  0:  BY  “ 


.  6  /  1 E6: AY  -  22 .  1 


Y  .  03 !  :  YC  -■  ..  1  1 3 
.6:  IF  MT  -r  70  THEN  GOTO  100 

o.zizr?.  —  ,  34:  E3  =  2.3:1-10  -  .000121  :X  =  i.4:XC 
053:3  =-  .  034:  AX  ~  -  4.0  /  1E6:  AY  =  79  /  1E6:BX 


1  ;£,f:  ;;  px 


.  61  : 
-i :  t:  . 


•  >  i v . vr 


:  I.r  MT  -  80  THEM  GOTO  1 OO 
<fl>  l:*  r  65:t:.Y  =  69:  PX  .  3:  ES  -  26.5:H0  ■=  .000125:  X  =  .  4:  XC  ~  . 

■4:5  ■-  .  .2 :  A X  -•  72.5  ’  IE 6:  AY  AX:PX  ~  On  BY  ~  A ■  jf  a 

GOTO  10O 


BY 


4  :  '  i~ 

I  IE  M 


1 

i  )t  > 

MR  -  1 

/  ■:  i  -  i 

X  *  PX  *  EY  /  EX) 

1 

1  0 

DIM  0 

<  3 ,  > 

1 

2>  • 

A  *  » 

IE.  tt# 

" 

1 

-  ( - 

o ;  i ,  i  ■> 

»  MF'  *  L 

X:Q(2,2>  -  MR  *  EX:  0(2 

,  1> 

= 

MF'  * 

1 

L 

4  1  1 

n  <  i  "  2  > 

=  0(2,1> 

5  0(3,3)  =  ES 

1. 

5  0 

d  -  (,!•: 

1  .,  1  ,■  i  0  : 

1 

6  <  j 

Ut  ■  ■■■■  ( 

3  *  0  +  2 

V  0 ' 1 , 2 )  +  4  *  c  (  3 .,  ( ) 

8 

1 

70 

02  ( 

D.  (  1.  .,  1  >  - 

0  <  2 2 )  >  /  2 

1 

90 

U3  t 

i  -  *  o 

(  1 , 2 )  4  t  Q ( 3 ,3) )  / 

8 

1 

90 

04  -  ( 

D  +  6  *  C 

(  !  .,  2  4  V  0  (3,3 )  >  / 

G 

fii  t 

1  IT  -■ 

7  2  i  t 

(1,2)  •  4  %  0(3,3))  / 

'is 

'4 

lo 

!  !  II  1.1 1 

"NUMBER 

OF  PLIES,  N— " ; M 

> 

i  i  ii  'p  i  i 

"  /II  1 

7FI'  GORE  !  Ft 1  0!  NESS  I  !> 

'U 

”  - 1 

ir 

H  M 

t  H"  /  •:  i 

HC)  :  HM  ==  (1  -  HP)  / 

N 

4  '0 

IF  HC 

0  THU'- 

M  =  N:  GOTO  260 

2 

5u 

INPUT 

"NUMBER 

OF  PLIES  BETWEEN  3  =  1./ 

2  AND 

CORE  , 

2! 

i  t 

PP I  M  T 

:  PRINT 

" PLY  ORIENTATION  (FROM 

2  zz 

FI. 

'2  TO 

7,  . 

DIM  0 

( 6  0  >  ,  P  ( 6  C 

,  4.  , X  (3,3, 9) , V  (4, 2)  .  W  1 

\  60 

l’1=  "  ;  M 
1  H  /  2  j  ’ 


J 
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.xth 

q,-> 


1.'  ~  T>I:  T 

F  DR  i  : 


~  1 
i  TO  N 


ZO'J 

IMPU 

i 

"p 

LY  ANCLE  = 

"  ;  0  ( 1 ) 

z  t 

IF  !M 

= 

IT 

THEN 

GOTO  330 

»;") 

IF  1 

IT 

THEN 

PRINT  "CORE, 

II 

;  I  NT  ( 

( 

--  N 

;  "  PLY 

THICK" 

33*’> 

0  ( I  > 

=- 

0  < 

I  >  * 

3  ..  1  4  i 

5926535  / 

1 30 

340 

P  [  1 , 0 

) 

= 

1 

350 

P  (1,1 

) 

- 

COS 

<2  * 

D  ( I  >  ) 

3  \j  \) 

P  ( 1 , 2 

) 

- 

cos 

(4  * 

0(1)) 

“* »“) 

P(I,  3 

> 

= 

SIN 

<  2  * 

0(D) 

330 

P  (1,4) 

as 

S I IM 

<4  * 

0(D) 

390 

IF  I 

M 

THEN 

D  = 

1 

400 

W 1  U  ) 

~ 

C: 

+  D  * 

HO  +  (  I 

- 

1  )  * 

HN 

4  i  0 

NEXT 

I 

420 

FOR 

K 

0  I'D 

2 :  K 1 

=  K  *■  1  : 

IF 

K  = 

2 

THEN  F 

430 

F  OR 

I 

=:= 

1  TO 

N 

<140 

IF  1 

~ 

0 

THEM 

TT  - 

HN 

450 

t  F  F 

- 

1 

T  HEN 

TT  = 

HN  *  (HN 

-t 

2  * 

W1 

( I)  ) 

460 

IF  ! 

= 

THEN 

TT  - 

(Wl(I)  + 

HN ) 

- 

W1  (  I 

470 

FOP 

J 

- 

0  TO 

4 

480 

3'  (  0  ,  F 

) 

= 

V  ( J  ,  I 

)  +.  p- 

(DO)  *  TT 

*  FF 

490 

NEXT 

J 

500 

NEXT 

I 

5 1  0 

Cl  '■= 

V  ( 

o. 

K  >  * 

Ul  :  C2 

=■=  9(1,  K  > 

* 

U2 : 

C3 

= 

9(2, 

2 ;  C5 

», 

( 4  ,  K  ) 

*  03 

;  .  2  *  /> 

■  n,  j 

,  F 

) 

=  Cl 

+  C2 

+■  03 

530 

Y  (  ]  .  2 

,  F 

> 

---  V  ( 0 

,  K )  * 

U4  -  C3 

340 

i  •:  l ,  7 

„  j. 

> 

=  C  4 

h-  C5 

550 

X  (  2  ;l  1 

,F 

) 

x  ( 1 

'  ">  l  -  \ 
<|  .  * 

560 

X  <  2  „  2 

,  F 

) 

=  Cl 

-  C2 

)  C3 

570 

v  ( 2 ,  3 

,  F 

) 

=■=  04 

-  C5 

580 

V  (  3  ,  | 

F 

) 

=  X  v  1 

,  3 ,  K ) 

590 

X  (  3  7 

•.  K 

> 

=  X  <  2 

,  3 ,  K ) 

600 

X  <3, 7 

,  l=: 

'I 

=  9(C 

,  i< )  * 

U5  -  C3 

6  1 0 

NEXT 

F 

*!!j  2  0 

L  I  ~ 

0: 

LO 

=■=  9: 

GOSUB  2310 

630 

F  =-  1 

;  LA 

»  9:  LB  =•  I 

» L.C  =•  7: 

GOSUB 

2420 

640 

F  »  1 

:  l. 

A 

»  IsL 

B  ~  9 

sLC  ■  8: 

GOSUB 

242 

0 

£j O 

F  - 

- 

3  2 

LM  — 

1. :  LB 

==  7 :  L.C  “ 

6: 

GOSUB 

'7> 

420 

660 

LA  - 

7 : 

LB 

6 : 

L..C  ^ 

6 :  GOSUB 

25 

20 

670 

!...  I  = 

£i  “ 

L.C 

6: 

GOSUB  2310 

r  = 

3  - 

LA 

7:  LB 

=•  6:  LC  « 

4: 

GOSUB 

71 

420 

67'.' 

F 

1  : 

L.A  = 

4:  LB 

=  S: LC  = 

o : 

GOSUB 

420 

7  '">!.> 

LA  =• 

O  • 

LB 

rr  T  • 

LC  ~ 

3:  GOSUB 

25 

20 

7 1 0 

FOR 

i 

r: 

1  TO 

’3 

j  2 » j 

FOR 

j 

” 

i  re 

7  30 

X  ( J  ,  I 

r: 

) 

-  X  <  I 

,0,4) 

7  40 

NEXT 

J 

750 

NEX  T 

T 

760 

IF  CN< 

:l ) 

-  o 

THEN 

GOTO  300 

M  (1  /  HN)  *  100  +  .  5 )  / 
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D*  "  : 


y-  .  i 

A  T 

-  "  B#  "  : 

V*  «  "  (G! 

■So 

X  1 

-  1  :  X2  -  1 

: !  ~  0 : 

GOSUB  2: 

XL.iM 

.<  1 

-  3 :  X  2  ~  1 

:  t  =  1  : 

GU63UB  21 

» ■  i 

!  f 

CM (2)  =  0 

THEN 

GOTO  040 

l  i> 

v  t 

"ALPHA#" 

:  Y 1.  =  " 

BETA*/3" 

(GPA) " : 1 T 


JB*  "  :  NT 


( 1/GPA)  “  :  :%■ 


LjGGUB 


i  RDETA*  *' :  NT  -  DEL 


:  GOSUB  2150 
82"  L  -  1  s  X2  *  i  /  T. :  i 
030  A  t  -  I  :  XI'  ~  1:1  =  5 


^  3:  GOSUB  . 
GOSUB  2220 


84'- 


i  i  i  !'  i  ■- . 


r>  111?  r !  GOTO  0630 


635 

XT  ■= 

"A 

"  » 

VT  - 

"  3  " :: 

,'T  -  »  "  :  I T 

r= 

"  B  "  :  WT  “  "  0  "  :  Gf 

'SUB 

'-i  ;  t:  ,*■ 

06'"' 

XI  - 

1: 

X2 

H 

#  H 

7  2:1,  =  0: 

GOSUB  2220 

S'?0 

X 1  = 

%  H 

?:  x 

2  -  H  3  / 

12 

:  t  - 

1 :  GOSUB  2220 

BOO 

if  l: 

N  ( 

4) 

■»  0 

THEN 

GOTO  920 

090 

XT  - 

"ALP 

HA"  : 

V  $  •“ 

"BETA":  XT  =• 

"TRBETA 

"  :  1*1  T  =■  "DEL..’ 

f'A"  : 

GOSUB 

v  o  o 

y  i 

1 

H:  X2 

<■1 

/  H  2:K  - 

3  s 

GOStUB  2220 

9  1  0 

X 1  ^ 

+  ’ 

H  ' 

2  s  X2 

-  12  /  H 

5  s  i-:-. 

5 : 

GQSIJB  2220 

#?  2;*> 

XT  = 

"  S3 

TR 

AIN 

AN  Ail  \ 

G IS” : YT  =  » 

"  : 

7  r  - 

"  " :  WT  “  ' 

" :  G 

OS  LIB 

930 

E  I M 

N  •' 

3 

1  >  ,  l-l 

",  1  \ 

, E (to, 12, 1 ) 

,st 

■•3,1) 

, 32 (3, 1 > , ET 

:  3  > 

940 

n  ( i , <: 

> 

1 

1£9;  |\| 

1 2 , 0)  ~  0 :  N 

<  3 

0  '•  ~ 

o 

950 

M  (  1  ,  C 

) 

:  M  ( 

2 ..  0 ) 

“  0 :  M  ( 3 , 0 ) 

=:  0 

V60 

F  OR 

I 

= 

1  TO 

9  70 

N  (  I  ,  C 

) 

- 

N  <  I  , 

0 )  / 

H 

930 

m  >:  i ,  v 

> 

=: 

M  ( I , 

0 )  * 

<H  #  H  /  6) 

990 

NEXT 

I 

l  ooo 

i  NF 

LJT 

II 

1  ENP 

.  DIFFERENCE  (IN 

K ) 

,  D  T 

";DT:  INF  LJT 

"  Mu  1 

GTURE 

_  II  „ 

PF 

INI 

1 0 1 0 

ET  ( 1 

) 

= 

AX  # 

D1  '■ 

BX  *  C : E  7  ( 

2 ) 

?=  AV 

*  DT  +  BY  * 

C 

1 02-: 

JJ  = 

1 

s 

IF 

ABS  ( 

DT  >  h  ABS 

(C) 

-  0 

THEN  JJ  =  0 

.1  O  3  C 

[  F 

J  J 

“ 

0  THEN 

GOTO  1130 

1  '040 

3.X  - 

G 

( l 

1  ) 

*  F  T  ■. 

1 )  +  Oil , 2 ) 

# 

ET  (2) 

J. 05  C 

SY  « 

Q 

\ 

,  1  > 

*  t:i  < 

1 )  o  0(2, 2 ) 

*: 

EZT  (2) 

1  OoL 

PM  - 

b 

#  (f, 

X  '  B 

t )  :  ON  -  ..  5 

#  ( 

SX  - 

SY) 

i  o  b: 

N  <  1 

1  ' 

PN 
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,0)  +  ON  # 

V  (  1 
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1  ) 
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V  ( 3 ,  1  > 

J  1.3C 

FT:  ■ 

m  r 

TAB  t 

a  >  ;  " 

EFFECTIVE  S 

TRE 

SSES" 

•1  * 

•EFI-E 

C 1  I VE 

1  i 4  L 

E  J  T 

0 ) 

■r 

"ME 

CHAM  I 

CAL  "  :  E.JT  (  1  > 

~ 

"NON- 

MECHANICAL" 

1  150 

A  $  \  1 , 0  > 

=  "N# 

AT (2,0)  = 

"M#  " : AT  (3,0)  =  " 

N 

1  160 

AT  (1,1) 

~  "MN*": 

AT (2,1)  * 

"MN#"  :  AT  (3,  1  >  ==  " 

NN 

1170 

F  OR  J  - 

0  10  JJ 

:  X  1  =■-  1 :  X  2 

'  a  i 

1  100 

PR  I  NT 

TAB (  2  ’  - 

; EJT ( J) 

1  190 

FOR  L 

0  TO  1  : 

PRINT  AT ( 

1  +  2  *  L  ,  J )  ;; 

1  200 

FOR  I  ■- 

t  TO  3: 

?  PRINT 

IJSEAT;  X  1  #  NCI,  J  ) 

? 

1  2  1.  0 

HE  X 1  I: 

PRINT  " 

" ; AT  (2 

*  2  #  L, J > ; : X j  a 

X  J 

1220 

FOR  I  -= 

l  TO  3: 

K.  PRINT 

LISE  AT:  X2  #  M  ( I ,  J  ) 

5 

1  2  30 

NEXT  I: 

X  2  =  X  2 

i  H  *  H  / 

6:  PRINT 

1240 

NEXT  L 

"N  " : A* (4, O'  =  "M 

A 4*  C  4 ,  1  )  =  "MN  " 


1 


1 

.  ‘i  >’ "• 

EOF:  P  = 
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.7  J 

1 

.V  o 

GOSUB  2 
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1 
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FOR  K  = 
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J 
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t 
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J  =  S 1  (I ,  J  > 

1 

E  O',  I 

3  *  L 

-i  6 ,  J  )  =  E  ( K 

t 

40 

ME  X  T  I 

1 

350 

NEXT  J 

t 

“T  i , ,"» 

FOR  J  - 

0  TO 

JJ 

1 

GOSUB  2 

660 

1 

ME  X  T  J 

l 

390 

NEXT  L 

1 

NEXT  F. 

1 

•!  i  o 

PRINT  : 

PRINT 

"PLY";  TAB! 

E  >  s  "LOWER  PLY  SURFACE"  ;  .  ;  "UPPER  »  ' 

L.Y  SURFACE" 

142-'  1!  C N  < 5 +  CM  -6)  -  0  THEN  GOTO  1560 

1430  71(0,0)  *•  "ON  AX  [S  MECHANICAL  STRAIN" : Z* (O,  1)  =  "OFF  AXIS  MECHANICAL 
STRAIN"  :  ;•:*  a  ,  0>  »  "ON  AXIS  NON -MECHANICAL  STRAIN"  : Z* ( 1 ,  i  )  *  "OFF  AXIS 
NON- MECHANI  CAL  STRA  [  N" 

iB  (  20  >  5  Z*  <  J  M6  >  :  GE  1  t.:' 


.!.  ‘VO 

EOF;  J  - 

■■>  7  0 

JJ 

145  ■ 

PF  1  NT  : 

FOR  1-1 

') 

-  1  ■■■  CN  ( 6 ) 

TO  CN(6>:  PRINT 

i  '!  ,<"• 

t  >  i  •••- 

}.  J  o 

1  i  •• 

print  I:  r. 

AB  (  5  > 

1  4 

F  OR  L  =■• 

o  to 

1 

j.  4  SO 

FOR  I  - 

1  TO 

3 

1  '■!  9 0 

PRINT 

USE AS ; 

E  ( K ,  I  +•  3  * 

L  +  M6  *  6 , J )  ; 

1  SCO 

ME  XT  1: 

PRINT 

•i 

5 

!.  5  1  ■ » 

ML  XT  L 

1520 

PR  I  NT 

J  530 

MI  XT  K 

1  540 
t  ‘.orv 

NEXT  M6 
ME  X  T  J 

i56o 

If  CN-.7) 

-  CM 

■,8 

>  -  0  THEN 

GOTO  2140 

1570 

D I  M  F  <  2 , 

2 ) ,  g  •: 

"T  ( 

3)  ,  EX (6) ,U( 

3)  ,  CM  (60.  1  )  ,  R  .1  (60 

1.5  SO 

1"  X  V  IS  T  AR  - 

-  •  > 

M  t.J 

.1  5  ?'"> 

F  <  1,1)  - 

1  ( 

L 

*  XC) : F (2, 2 

>  =  1  /  (Y  *  YC) 

1 6 

F-  (  1 , 2 )  - 

F  XYST 

;TP 

*  8 OR  (F ( 

1,1)  *  F(2,2> ) :F< 

1  6  1  0 

FX  -  1  / 

X  1 

/ 

XC: I Y  =  1 

/  Y  -  1  /  YC 

1  620 

EX (1>  -  X 

/  EX 

:  E 

v  ( '2 )  =  y  / 

EY; EX (3)  =  S  /  ES 

-  Y  C 

/  EY: 

EX 

(6 >  ~  -  S 

/  ES 

1  630 

py  =  FX  T 

0(1, 

1 ) 

+  FY  *  0(1 

,  2) 

16  40 

GY  "  FX  * 

0(1, 

2) 

+  FY  *  0(2 

,2) 

1650 

FOR  K  - 

1  TO 

1 660 

FOR  F  = 

1  TO 

!  670 

0  ( 1  ,  F  )  = 

(j 

160" 

FOR  I  = 

1  TO 

1690 

FOR  J  - 

1  TO 

!  700 

bir  ,F!  = 

G  ( K ,  F 

\ 

+  F  (  I  ,  J  >  # 

0  <  I ,  K )  *  O  ( J ,  F ) 

1  7 1  0 

NEXT  J 

17  20 

NEXT  I 

1730 

NEXT  F 

1740 

NEXT  I 

1  750 

G  (  7,3/  “= 

(0(3, 

'i 

/  S)  2 

R  ( 5 ) 


F  (  1 


XC 


EX::  EX  (5 


18 
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1/6';  FOR  :  -  !  1  i  !  K 
1  .•  I  t  if-:  L  -  "  ro  1 

l  7 30  A  ~  0:H  -=  O :  C  - 


;  '90 

F  OR  1 

-  1  1  U  ’■ 

:  FI  I  43*  L:BB  = 

0 

1  300 

1  e  •:  1 

,  I L  ,  0  > 

O  THEN  BB  -  3 

1  E.<  1  ■  ? 

J  <  1  :  --- 

4..  1  1  ’ 

7  •:  f  ,11.  -,  1 

ET  (  1  )  :  IF  Eit  ,11. 

,  0 )  “  0 

THEN  ECi,  1L,  Ol  =  F  (  F.  „  I 

,  320 

R  l  I  >  =■ 

1  p  >,  <  T  1 

£i.O  >  -  U  ■:  J  )  >  /  Ft!.,  i  L 

, 0) :  IF 

I  ’  .1  rHETl  GOTO  iU4o 

1  R30 

[  r  r  <  1 

)  RNO 

,  L)  Ti  iEIM  GOTO  1050 

iv'OSO 

RM  ( 1-  ,  L) 

----  R  •  I  )  ; 

'.Mil  ,  L)  =  I 

1 950 

NEXT  r 

1  2  tj  0 

FOR  I 

-  1  TO  7 

:  .1. 1  -  1  +  3  *  L 

1070 

FOP  J 

-  1  TO  7 

:JL  =•  J  4  3  *  L 

1  BOO 

A  -  A  ' 

G  C I ,  J  > 

*  E£(K,IL,0)  *  ECi.  JL. 

,0):  IP 

JJ  -  0  THEN  GOTO  1910 

1  390 

u  =■  t:.i  + 

G (I , J ) 

*  CE  <K„  IL,0)  *  LI  <  J  ) 

4  E  C  K  ,  J  I. 

,  0 )  t.  LI  C  3  )  ) 

J  900 

c  *  c  + 

G (I , J > 

*  U  (  I  )  *  U  v  J  ) 

1  9  i  0 

NEXT’  J 

1 920 

NE’X  1  i 

1  93,.. 

D  "  B  4 

GX  *  E  C 

1,1  +  3  *  L. ,  0)  +  GY 

*  E  <  K ,  2 

4  3  4  L ,  0 ) 

1  9  40 

c  ~  C  4 

OX  #  L-  C 

1 )  +  GY  *  L)  C  2 )  -  1 

.!,  950 

V  2  S 

QR  ( B 

2  4  *  A  *  C) 

1  96*0 

R  i  C  H  ,  L  ) 

-=  (  -  [• 

4  02)  /  2  /  A 

1  770 

NEXT  1... 

1  90'..' 

NEXT  1. 

1  99,, 

X  1  --  "  s 

150 

I'RENGTH 

ANALYSIS" :Y*  »  ’’  ":V 

f  ~  n  H . 

Z*  -  "  "  ;W$  =■  "  GOG 

- 

FOR  I 

•=  1  TO  2 

2<.»  i  0 

PRINT 

II  It  .  II 

TSAI -HILL"; "  ";"MAX 

.  SIR."; 

"  " ; "ST R .  COMP."; 

2020 

NEXT  I 

2  o  3  0 

F  E 1  NT 

::  PRINT 

20  a  0 

X  1  =■-■  H: 

Nil)  - 

"STRENGTH  RAT  10,  R"; 

’  i  c 0  >  » 

"NORMALIZED  STRENGTH  !■ 

,  R .'  H  ” 

2u 

r  or.  k 

-=  1  O 

(7)  TO  CN<8> 

'' 

PRINT 

TAB  C  20 

)  ;  FTt  NO  :  PRINT 

*  70 

E;  OR  1 

■=  1  I'O  N 

PRINT  I ,  TAB  v  5); 

OR  L  *  O  TO  J 


ZO70  9  --PINT  USEASsRl C I, L>  /  XI;:  PRINT  "  "  ;  ;  &  PRINT  US  EH!  At ;;  R!i  <  I ,  1 

XI;  :  PR  INC  "  "  ;  CN  ( I ,  L)  ;  ;  .  ■ 

next  l 

2 1 10  PRINT 
2120  NEX  T  1 

2 1 30  •>.  1  -  L :  GE!  G*:  NEXT  K 
.■1  4','  PNC) 

215''  PRINT 

*  ’  !  t.o  PR  I  NT  "******  i *  *  4  *  *  *  *  *  *#  #  *  *  !f  ********  *  *  *  *  *  *  *  *  *  " 

:.'!/>•  PRINT  TAB  <  15)  ;X*;"  Y*;  "  V$ 

2  IQ’.'  PRINT  TAB  v  15)  ;  Z*;  ”  " ;  W* 

'2 IV"  G<:  !  Di 
2200  PRINT 
2.2  i  o  RETURN 


19 


1 


r 


. .  ,  1  . 

.  “1.  i  i 

TO  3: XM  = 

X  1 

-  2  **>’,> 

!•  OR  L  -  t 

TO  h  *■  I : 

IF  L 

2.  1.'  4  *  1 

TOR  J  -  i 

TO  3 

-  .  '1  C"  , 

PRINT 

USEAt; XM  * 

x  a. 

=  K  +  1  7  HEN  XI'I 
,  L)  : 


:  '.'f-  next  u 
227 V  NEXT  I. 

"'230  PRINT 
229 NEXT  I 
2 3 VO  RETURN 
. 3  1 0  ft  =  X  •.  1 ,  1 ,  L I  >  :  B 
•  F  X  <  3 3 ,  L.  I  ■ 


X  <1, 2 , L  I 


U, 


.  I )  :  D 


x  t : 


*7  .  •, 

DF 

T  =  ft  * 

iD  *  F 

-  E 

* 

E) 

+ 

B 

*  (2  *  E  *  C  - 

F  *  B»  -  D 

^ .  .j  S  ’  . 

i 

1,1, LO) 

*  a.)  * 

C 

c 

* 

E) 

/ 

DET 

::%34  o 

2  "'St./ 

Y 

V 

1 , 2, LO) 
.2,  1  ,  LO) 

<0  $ 
-  x  >:  i , 

E  ~ 

2,  L.0) 

B 

* 

F  > 

/ 

PET 

.. :  cJ 

X 

x 

I , 3 , LG ) 
23,  1  ,  LO) 

~  i  B  # 

-  XU, 

E  ■ 

3 ,  LO  > 

D 

* 

C) 

/ 

DET 

2380 

v 

2,  L0> 

=  ■:  a  * 

F 

F 

* 

C> 

/ 

DET 

"T  ) 

x 

3,ld> 

-  <  B  T 

C  - 

A 

* 

E> 

/ 

DET: X (3, 2, LO) 

-  X  (2, 3,1.0) 

"•i\  or- 

24  1 0 

r 

3 ,  3 ,  L  Q  > 
UPTURN 

.r:  A  * 

D  -- 

B 

* 

D) 

/ 

DET 

*  c  *  c 


••J 


■or.  i 

■OR  J 


TO 

T  0 


- '  4  4  0 

ru  ■=•  o 

245'./ 

FOR  K  * 

l  TO  3 

2460 

6LJ  -  SU  + 

X  ( I  .  K 

,  LA 

) 

% 

X  (K,  J,L 

D> 

2  4  '0 

NEXT  r 

24  E" 

X(I,J, LC ) 

=  F  * 

SU 

2490 

NEXT  J 

2:5  oO 

NEXT  I 

25 1  0 

RETURN 

25 

-■■OR  I  = 

1  TO  3 

2530 

FOR  ,J  « 

1  TO  3 

25  4  0 

X  •:  i ,  j ,  L.c: ) 

=  x  ;  i 

,  J. 

LA  > 

+  X  ( I  ,  J 

,  LB) 

2550 

NEXT  J 

T  6" 

NEXT  I 

2' 5  .  20 

RETURN 

;■  i 

FOR  I  ■* 

I  TO  3 

PTi/PO 

5  KI,  P  >  -= 

0:  S2  ( 

I ,  v- 

) 

= 

o 

260') 

FOR  J  = 

1  TO  3 

2'  i  i  0 

5 1  •:  i ,  p )  =■ 

Sl(l, 

P ) 

+ 

X 

(I,  J,3> 

* 

■7 

£ 

~T1 

>  o ...  ' .  ’ 

32(1,  P>  =» 

32(1, 

P) 

+■ 

x 

( I ,  J  ,  5  / 

*  N  ( J  ,  P ) 

26  3/ 

NEXT  J 

2  6  40 

NEXT  1 

26  54  ■ 

RETURN 

.’6  60 

T'  ~  .5  * 

(E  U  .  1 

+ 

* 

L ,  J  >  + 

E  ( K ,  2  + 

26  70 

0  -  .  5  * 

(E4  ,  ! 

-i 

T, 

* 

r 

C-j 

! 

E  ( K ,  2  + 

A  V  L 

X  (  I 


*  L, 

*  L  , 


,  )  .4  \ 

M  ‘.J  if  .  f 

,  J  ,  6  ) 


i )  > 

J )  ) 


m  •:  J p ) 
M  ( J  P ) 


263'/ 

2690 


!': 

E  ( 

.5  *  E  < K ,  3  +  3  * 

K ,  1  +■  3  *  L. ,  J  >  «  P 

L, 

+ 

J) 

Q  * 

P  ( K ,  1  )  + 

R  * 

P  ( K  ,  3  ) 

E  1 

'-:,2  +  3  *  l,j>  =  r> 

- 

Q  * 

P(K,i>  - 

R  * 

P  d  , 

1,3  +  3  *  L.,4)  2 

* 

R  * 

P  ( k ,  1  )  - 

2  * 

u  t  <  k: ,,  .3 ) 

RF  T I  iRt  i 
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46E  ‘  vm.) 

r),  J  i  OB 

l 

2 .  9  ■  ' 

J  ..  V  7 1  t 

'Ll.  +  on 

:  ..  9‘jE.  1  09 

1  .  V  PE:  *  08 

.'  .  7  A  * 

v  .i  t: '  of: 

0.2/9  '9 

yt 

1 . 3  Mil  i  .0, 
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1 , 9d»E  *■ 
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1  .  . 

i  e+  on 

2. 02P  1  ■  9 

T 

2 ..  1  1  I  :'  i  •  ±; 

!  ..  UAL  1 

OOE:  -OB 

:  .  a  O'  E  V  .’S 

1 . 92L  +  '.'9 

1  ..  6 ML  J- 

L TRENGT H 

RA T ! U , 

L: 

7  :•  r:  +  os 

t  .  n9L 

1 

L  28! 

■1  ..  4  6 b  • 

■  ■  0  ■  .  L 

,_,Cy[  ..n'.j 

1 

1  .  '  -  71  i 

I  99L> 

l  HE:  *  ’ 

.1:  .sil  t-  05 

2. o.'.r-  ••.t 

' ..  2  l  L  4 

1  .97  ’  •  ’H 
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Problem  #3 


AMU  . .  AL  I  -T-’f:,  '!  ~ '3 1 ' 

mum:  !  ur  i-T.  r i.  s ,  im 

HI ;!  .''i  H.  I  ZED  CORE  THICKNESS,  HLVH=U 

r  .  OR I ENT AT  TOM  >.  FROM  l~-  H/J  TO  Z  ;'*H  '  M  > 
PL  ,  ANGLE  “O 
r  L  .-  ANGLE  O 
Mr  ANGLE  “O 
Pi..  /  ANGLE  "O 


!  a  mule: 

:  i  angle:  o 

F  :  v  ANGLE  *  ,!C 
Li  .  ANGLE:  =90 
!■  i  /  ANGLE  90 
1  -LV  ANGLE  =90 
PL/  AUG ME  90 
ANGLE.  =90 
PLY  ANGLE  90 

,:T.  Y  angle:  9 <' ■ 

?  *  *  4  r  *  *  %  %  *  *  \  #  *  *  t  t  *  *  *  *  t  ♦  %  *  *  *  *  #  I  ■*:  1  >  *  *  *  4 
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